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Within the framework of a modified proton ordering model of the KH2PO4 family ferroelectric crystals, taking 
into account a linear over the strain £6 contribution into the proton system energy, we obtain an expression 
for longitudinal dynamic dielectric permittivity of a mechanically clamped crystal using the four-particle cluster 
approximation and the dynamic Glauber approach. At a proper choice of the model parameters, we obtain a 
good quantitative description of available experimental data for these crystals. 
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1. Introduction 

In the late 1960-ies, most theoretical and experimental studies of ferroelectrics concentrated on var- 
ious dynamic phenomena. The dispersion of dielectric permittivity of ferroelectrics was explored at low 
frequencies, which provided an important information on the mechanisms of phase transitions and re- 
vealed the peculiarities of the low-frequency dynamics of a system. Ferroelectric dispersion is closely 
related to the presence of a low-frequency excitation, i.e., a soft mode which can be either resonant or 
relaxational. Ferroelectric compounds of the KH2PO4 family occupy an intermediate position. The region 
of fundamental dispersion in these crystals is located in the submillimeter range v « 50 GHz. At deuter- 
ation, the ferroelectric dispersion in these crystals is shifted to the millimeter and microwave ranges. 

The major task of dielectric spectra studies of ferroelectric crystals is to explore the peculiarities of the 
soft mode behavior, especially in the phase transition region (lfl. As a rule, the soft modes in the KH2PO4 
family ferroelectrics are strongly damped. To explore their character is a complicated task. One has to 
explore the dielectric spectra of these crystals in a wide frequency range that includes several regions re- 
quiring specific and unique experimental methods of measurements. There is hardly any experimental 
group fully equipped for such studies. This fact, along with the principal difficulties in experimental mea- 
surements of dielectric spectra, and the dependence of e* {u>, T) on sample quality and surface treatment, 
causes the situation when the experimental data for dielectric spectra of the KH2PO4 family ferroelectrics 
turn out to be disembodied and quite conflicting. This should be kept in mind while analysing the exper- 
imental data and the theoretical results for dynamic characteristics of ferroelectrics including those of 
the KH2PO4 family. 

In the late 1970-ies, the obtained experimental results for the dynamic characteristics in the KH2PO4 
family compounds were interpreted mostly within phenomenological models (see Hi). Phenomeno- 
logical theories do not make it possible to reveal the microscopic nature of the dispersion of dielectric 
permittivity or to appropriately describe the effect of various factors on the character of its temperature 
and frequency dependencies. The attempts to solve this problem using the Green's function method or 
Bloch kinetic equations method failed 
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Avast majority of studies on the theory of relaxation phenomena in the KH2PO4 family ferroelectrics 
are based on the stochastic Glauber model |6]. For the first time, the relaxation dynamics of the KD2PO4 
type ferroelectrics was studied using this method in £1, where, within the four-particle cluster approxi- 
mation (FPCA), there was initiated a study of the main regularities of longitudinal relaxation in the case 
of a paraelectric phase. However, long-range interactions were not taken into account therein, and the 
corresponding experimental data for the KD2PO4 type ferroelectrics were not discussed. Later on 
a more consistent model of deuterated KD2PO4 type ferroelectrics and ND4D2PO4 type antiferroelectrics 
was explored. Within the framework of this model, using the FPCA for short-range interactions and the 
mean field approximation for long-ran ge in teractions, longitudinal dynamic characteristics of these crys- 
tals were calculated. It was shown I10l4l2ll that the theory proposed in provides a satisfactory de- 
scription of thermodynamic and longitudinal dynamic characteristics of the KH2PO4 type ferroelectrics. 
In fl3Hl5ll . the authors attempted to develop a more consistent theory of the KH2PO4 family ferroelectrics 
in the FPCA which takes tunneling (Q) into account. The results were not good enough to appropriately 
describe the available experimental data for the dynamic characteristics of these crystals. However, the 
fact of suppression of the dynamic characteristics of the KH2PO4 type ferroelectrics by short-range inter- 
actions was established. An effective tunneling parameter Q (O « O) renormalized by the short-range 
interactions was obtained. It should be noted that the established in flilJisll suppression of dynamic char- 
acteristics of the KH2PO4 type ferroelectrics by short-range correlations is the most probable reason of 
the De bye-typ e dispersion of dielectric permittivity observed in these crystals. 

In I16l4l8ll . thermodynamic and dynamic characteristics of quasi-one-dimensional hydrogen bonded 
CSH2PO4 ferroelectrics were found using a self-consistent approach to the calculation of thermodynamic 
and dynamic characteristics of pseudospin systems with essential short-range and long-range interac- 
tions, based on the calculation of the free energy functional with short-range interactions taken into 
account in the reference approach. It was established that an essential suppression of the soft vibration 
mode by short-range correlations takes place in a wide temperature range. This fact, just like in the case 
of KH2PO4, is directly related to the Debye type of longitudinal dielectric permittivity dispersion observed 
in CSH2PO4. It should be mentioned that similar studies of thermodynamic and dynamic characteristics 
of KH2PO4 can be carried out using the technique developed in 11911 . Such studies would make it possible 
to explore the effect of suppression of the soft mode in the KH2PO4 type ferroelectrics more consistently 
than in Hill and thereby to explain the Debye character of the dielectric permittivity dispersion in these 
crystals. 

It should be noted that the ferroelectric compounds of the KH2PO4 family are piezoelectric. Piezoelec- 
tric coupling is observed in external electric fields and mechanical stresses of certain symmetries. Ferro- 
electric phase transition in the KH2PO4 type crystals is accompanied by the appearance of spontaneous 
strains, which changes their tetragonal symmetry. So far, the calculations of dielectric characteristics of 
these crystals within the proton ordering model were restricted to a static limit and high- 

frequency relaxation. The attempts to explore the piezoelectric resonance phenomenon within a model 
that does not take into account the piezoelectric coupling were vain. The conventional proton ordering 
model does not permit one to describe the effects associated with the differences of the free and clamped 
crystal regimes in the static limit or the phenomenon of crystal clamping by a high-frequency field. This 
leads, in particular, to some quantitative deviations from experiment for the temperature behavior of 
polarization relaxation time and dynamic dielectric permittivity of the KH2PO4 type ferroelectrics in the 
phase transition region. 

The studies of the piezoelectric coupling effect on the phase transition and on physical characteristics 
of the KH2PO4 type ferroelectrics were initiated in(2oll, where the Slater theory (illl was modified by 
taking into account the splitting of the lowest ferroelectric level due to the strain £6. 

The most fundamental results for the KH2PO4 family ferroelectrics were obtained in HBJioll . For the 
deformed crystals of the KH2PO4 type, the Hamiltonian of the proton ordering model was modified for 
the first time by including £6 into in the shear strain (2^)2^], taking into account the deformational mean 
field and the splitting of the lateral proton configurations. Later on (24} [2^], all possible splittings of pro- 
ton configurational energies by the strain £6 were included into the model. In 12411 , using this model, the 
phase transition, thermodynamic and longitudinal dielectric, piezoelectric, and elastic characteristics of 
K(Ho.i2Do.88)2P04, as well as the effect of the oe on these quantities were explored. The same character- 
istics for other K(Hi_ x D x )2P04 type ferroelectrics were later on calculated in E^l . 
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The thermodynamic and longitudinal dielectric, piezoelectric, and elastic characteristics of the 
KH2PO4 type were also calculated in (25], |27|] within a model that takes into account the tunneling and 
piezoelectric coupling. It should be mentioned, however, that taking into account the tunneling within 
the cluster approximation yields a non-physical behavior of the calculated quantities at low tempera- 
tures lilll . In Eii - tioll . the effect of the electric field £3 on the phase transition and on the physical char- 
acteristics of K(Ho.i2Do.88)2P04 and KH2PO4 was explored, and a good agreement with experiment was 
obtained. 

In E3 - l28ll . where there was used a model with tunneling, the dynamic characteristics of the KH2PO4 
family ferroelectrics were not considered. In liill . using a modified proton ordering model proposed 
in 12411 , the dynamic dielectric permittivity of a free KH2PO4 type crystals was calculated taking into 
account the dynamics of £6 strain. The experimentally observed effects of crystal clamping by a high- 
frequency electric field and piezoelectric resonance in KH2PO4 and KD2PO4 crystals were theoretically 
described for the first time. Peculiarities of the ultrasound attenuation coefficients near the phase tran- 
sition temperature in these crystals were also described. In 13311 , we presented a detailed review of the 
obtained results for longitudinal and transverse static dielectric permittivities, for piezoelectric coeffi- 
cients, and for elastic constants of several ferroelectric crystals of the KH2PO4 family. Moreover, the 
typical behavior of longitudinal and transverse characteristics of mechanically free KH2PO4, Rb2P04, 
KH2ASO4 crystals was shown and the results for temperature and frequency dependencies of longitu- 
dinal and transverse dielectric permittivities of KH2PO4 were presented, along with the corresponding 
experimental data. 

In the present paper, using the model proposed in jzill we calculate the longitudinal dynamic dielec- 
tric permittivity of clamped ferroelectrics of the KH2PO4 type and explore its behavior in wide temper- 
ature and frequency ranges. Using the obtained results, we perform a detailed analysis of the available 
experimental data for these crystals. 



2. Systems of equations for the time-dependent deuteron distribution 
functions 

We shall consider a system of deuterons moving on the O-D. . . bonds in deuterated KD2PO4 type fer- 
roelectrics. A primitive cell of the Bravais lattice of these crystals consists of two neighboring tetrahedra 
P0 4 along with four hydrogen bonds attached to one of them (the "A" type tetrahedron). The hydrogen 
bonds attached to the other tetrahedron ("B" type) belong to the four structural elements surrounding 
this tetrahedron (figure [1). 




Figure 1. A primitive cell of the KD2PO4 type crystal. One of the numerous possible ferroelectric deuteron 
configurations is shown. 

The dynamic characteristics of these compounds will be calculated within the four-particle cluster 
approximation that proved to be successful in describing their thermodynamic properties fl2lE4T - l27ll . 
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The Hamiltonian of the deuteron subsystem, taking into account short-range and long-range inter- 
actions in the presence of an external electric field £3 along the crystallographic c axis and mechanical 
stress 06 - <Jxy, which independently contribute to polarization £3 and strain £§, consists of the "seed" 
and pseudospin parts I24H2H1 : 



H = NH m + K 



(2.1) 



where N is the total number of primitive cells. The "seed" energy of a primitive cell corresponds to 
the sublattice of heavy ions and does not explicitly depend on the deuteron subsystem configuration. 
It is expressed in terms of the strain £6 and electric field £3 and includes the elastic, piezoelectric, and 
dielectric contributions 

1 „£0„2 " - 



H 



(0) 



£0 r-2 



(2.2) 



where v is the primitive cell volume; Cgg , e^g, X33 are the "seed" elastic constant, piezoelectric coefficient, 
and dielectric susceptibility, respectively. They determine the temperature behavior of the corresponding 
observable quantities at temperatures far from the phase transition T c . 
The pseudospin part of the Hamiltonian reads 



q'f 



£2^ 6 £ 6 ^- 



if z 



(2.3) 



The first term describes effective long-range interactions between deuterons; a q f is the z-th component 
of the pseudospin operator that describes the state of a deuteron in the g-th cell on the /-th bond. (/ = 
1,2,3,4). Two eigenvalues of the operator a q f - +1 correspond to two possible positions of the deuteron 
on the bond denoted by "1", "2" in figure[l] In (2.3) H S h S (.6) is a linear over the strain £6 Hamiltonian of 
the short-range interactions between deuterons jj^l : 



tfsh.s(6) = 



E 



8s6 $16 ), -> 

— £6- — £6 I (o>i + a q2 +G q3 +a q4 ) 

— £6 + ~ £ 6 I \P q\ a qlO qS + O qlO q2<J qA + V q\P q3<? qi + & q2<J qlO qi) 

+ j (V s + S a6 e 6 ) [<Jq\Oq2 + OqlOqi) + J (V s -<5 fl6 £ 6 ) (o q2 G q 3 + O qA (J q \) 
1 . ,1 

+ -U s [cr q ia q 3 + o q 2a q 4j + —®sO q \o q2 a q3 a qi 



(2.4) 



Here 



and 



1 1 

V s = —-Wi, U s = —Wi—e, <t> s -4e-8w + 2wi 



£ = £ a -£ s , W = Ei-E s , Wi=£ -£s 



where £ s , e a , E\, eq are the energies of deutron configurations near the PO4 group. 

The third term in (2.3) is a linear over the shear strain eq mean field Hamiltonian induced by the 
piezoelectric coupling; y/s is the parameter of the deformational mean field. 

The last term in (2.3) effectively describes the interactions of deuterons with an external electric field 
£3. Here fif$ is the effective dipole moment related to the f-ih hydrogen bonds, where 

^13 = P23 = A*33 = P43 = M3 . 



and ^3 is the dipole moment of up/down deuteron configurations. 

Taking into account the peculiarities of the crystalline structure of the MD2XO4 the type ferro- 
electrics, their dynamic characteristics can be calculated within the four-particle cluster approximation 
that proved to be effective in describing the thermodynamic characteristics of these crystals E2IE4II27I1 . 



33705-4 



Longitudinal relaxation of mechanically clamped crystals 



Long-range interactions are taken into account in the mean field approximation. Within the cluster ap- 
proach, the thermodynamic potential of MD2XO4 ferroelectrics calculated per one primitive cell reads 

g s m = H (0) +2v c (r^f + ^ £ InZ™ - ~ln^ - va 6 e 6 , (2.5) 

Z P f=i P 

where 4v c = /n(0) +2/12(0) + /i3(0), the eigenvalues of Fourier-transform of the long-range interaction 
matrix J ff /= £ Jff(qq'); 



R 9 _R r,' 



V (1)Z = (0>1> = <0-<?2> = <0>3> = (o>4> 



m -fiH^ (a) _/jr>( 4 ) 1 
is the parameter of deuteron ordering; - Sp e if s ,Zg s -Spe ^ <? s , p = are the single-particle 

and four-particle partition functions. The single-particle H^j s and four-particle H q ^ deuteron Hamilto- 
nians are presented by 

— ■ (2.6) 

lf s p 2 



2 



(2.7) 



where 



z 6 = jS (-AS + 2v c r/ 1)z - 2t// 6 £ 6 + /i 3 £ 3 ) , z 6 = /3 (-2A= + 2v c f7 (1)z - 2^ 6 £ 6 + ^£3) 



The effective field A£ exerted by the neighboring hydrogen bonds from outside the cluster, is determined 
from the self-consistency condition: the mean values (<J q f) calculated within the four-particle and one- 
particle cluster approximations should coincide. 

The dynamic characteristics of the MD2XO4 crystals will be explored usin g the p roposed dynamic 
model based on a stochastic Glauber model @], Using the method developed in llol - [l2l[32ll . the system of 
equations for the time-dependent deuteron distribution functions is obtained in the form 



r 



l-a qf ,tanh-pe z qf ,(t) 



(2.8) 



where £ q ^,(.t) is the local field acting on the f'-th deuteron in the q-th cell, which can be obtained from 
the Hamiltonian (2.3). Expanding tanhj/te^Cf) over the pseudospin operators <r q f, occurring in Hamil- 
tonian (2.3), taking into account the fact that o q f = +1 and the symmetry of the deuteron distribution 
functions in the MD2XO4 ferroelectrics in the presence of the electric field £3 



,d)z 



,(3)z 



= {Oq\) = {o q2 ) = {(T q3 ) = ((7, 4 ), 

= <0 q iO q2 - q3 ) = {(J q \(J cfiCJ qi ) = ((J q i(J q2 (J qi ) = {cr q2 cr q3 o qi ) 

= (o q2 o q3 ) - {a ql a q4 ) , tj 2 ={a q \(J q2 ) = {(j q3 (j qi ), 



((Tql(7q3) = ((rq2(rq4), (2.9) 



from (2.7), one can obtain a closed system of equations for the time-dependent single-particle, three- 
particle, and pair distribution functions of deuterons in MD2XO4 within the four-particle cluster approx- 
imation and for a single-particle distribution function within the single-particle approximation 13211 : 



dt 



^(3)z 

ft 

v if )z 



I Cn C U C13 C U C15 ^ 

C21 C 22 C 2 3 C 24 C 25 

C31 C 32 C33 C34 C35 

C41 C42 C43 C44 C45 

V C51 C52 C53 C54 C55 J 



^(3)z 
iff* 

v? z 
v vf )z 



I ci ^ 

c 2 
£3 
c 4 
\ c 5 J 



(2.10) 
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where the following notations are used 



Cll 


= -(i-P 6 z -Qli-Qffi) 




C12 = 


*6. c 13 = Ml v 


C14 = 


Kt- 5 i5 = Nl, 


ci = L^, 


C21 


= (2P| + 2Q| 1 + 2Q| 2 + 3i? 6 ) I 


C22 = 


-(3-^6 Z -Q61-Q62). 


C23 = 


UV 6 Z + M 6 Z 2 + L Z ), 




C 2 4 


= UV 6 Z + M Z 1+ L Z ), 




C25 = 


(M Z 1+ M 6 Z 2 +L Z ), 


C2 = 


UV 6 Z + M Z 1+ M Z 2 ), 




C31 


= 2UV 6 Z + M Z 2 + L Z ), 


C32 = 


2M 6 Z !, 


c 33 = -2(l-i? z ), 


C34 ; 


= 2P 6 Z , c 35 =2Qf 1 , 


c 3 = 2Q z 2 


C41 


= 2UV| + M ( f 1 +L z J, 


C42 = 


2M 6 Z 2 , 


c 43 = 2Pq, c 44 


= -2(1 


_ ^6)' C45 = 2Q| 2 , 


c 4 = 2(& 


C51 


= 2(M 6 Z 1 + M 6 Z 2 + L Z ), 


C52 = 


2JV£, 


£53 = 2(21!, C54 


= 2Q 6 Z 2 , 


C 55 = -2(1-.R Z ), 


c 5 = 2P§, 



(2.11) 



a— 77 (1)z = -?7 (1)z + tanhiz 6 . (2.12) 



3. Relaxational dynamics of mechanically clamped MD 2 X0 4 crystals 

Now, using the obtained systems of equations, let us calculate the dynamic characteristics of the 
MD2XO4 crystals. Let us consider the case of small deviations of the considered system from equilib- 
rium. We can separate the static and dynamic parts in the obtained system of equations. To do so, we 
present the distribution functions and the effective fields as sums of the equilibrium functions and their 
fluctuations 



,(3)z 



z 6 = z 6 + z 6 t, z 6 = -j8A£ + 2/3v c ti 



(i) 



(2)z 



■rjfWT tf = 1,2,3), 

zat = -PA c st + 2?>v4] )z + P^E 3t 



(3.1) 



Owing to a piezoelectric coupling, time-dependent electric fields should induce time-dependent strains. 
However, in the present paper we shall consider the fields with the frequencies of the order of 10 9 ~ 
10 12 Hz, which is far above the frequency of piezoelectric resonance. When the frequency is that high, 
the strains are not capable of following the external fields, which means that the crystal is effectively 
clamped. Therefore, in the expansions (3.1) we assume the strain £6 to be time-independent. 

We expand the expressions for the coefficients P z , . . . , Lg in series in Ze t l2 up to the linear terms. Tak- 
ing into account these expansions and (3.1), we obtain a system of equations that describes the behavior 
of fluctuational parts of distribution functions 1 3 2[ 1 3411 : 



_d 

dz 



( ^ 1 




( c n 
C21 


C12 
c 22 


Cl3 
C23 


C14 
C24 


Cl5 ' 
C25 


1 






C31 


C32 


C33 


C34 


C35 




'l2t 




c 4 i 


C42 


C43 


C44 


C45 




„(2)z 

v Vst > 




I c 5 i 


C52 


C53 


C54 


C55 ) 


V 



where the coefficients of the system read 

C24 = I(^- fc f^). 



C\2 ■ 



Cl5 = 



C22 ■ 



-J. 



C25 = 



c (0) 

L 12 




L 12 


c (0) 




c (0) 






c (0) 

u 12 


c l0) 

L 25 




c m 



H 3 E 3t 
2kT 



C\ 

C2 
C3 
C4 
C-5 



Cl3 



C23 



; (0) 
'13 



; (0) 

'23 



c 2 - —k 
a 



(3), 



K s C 13 



1.(3) =(0) 



(3.2) 
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cn = l(c™ + Pv c Y®- k Vv*), 

r _ 1 f~(0) 1.(2) -(0)") 

34 ~ a I 34 ~ sl 14 J' 

r _ 1 |=(0) 1.(2) =(0)1 
a V 44 ^2 C 14j- 

^ _ 1 |=(0) j.(2) =(0)1 
ffiV 54 53 14 J' 



C32 : 
C35 : 
C42 : 
C45 : 
C52 : 
C55 : 



=( 
=( 
=( 



SO) _ h {2) 
'32 



2) =(0)1 
1 L 12 J 



r (0)_j.(2) = (0)1 
L 35 Nl L 15 j 

■M _ 1.(2)^(0)1 
'42 ^s2 L 12 



-M _ 1.(2) =(0)1 
'45 ""52 L 15 J 

;(0) _ 1.(2) =(0)1 
'52 ^s3 u 12 



r (0)_ 7.(2) -(0)1 



'55 ,v s3 



a J 



C3 = -^V S , 

a 51 



(0) _ J-( 2 M°>1 



C43 ~ a I u 



C4 - — fc.o ?6, 



r (0)_ 1.(2) =(0)1 



C53 ~~ z ^ 5 
c 5 — — k[J r s . 



a 



The expressions for the quantities entering the coefficients of the system (3.2), are given in 13411 : if the 
piezoelectric coupling is neglected, they coincide with the corresponding expressions of Hill . 

The system of equations (3.2) is reduced to a non-uniform differential equation with constant coeffi- 
cients for a single-particle distribution function 



(l)z - (l)z -(l)z -(l)z ■ (l)z m . UsEzt „ 

It +P* T 't +P3lt +P2^t +Pl1 t +P0^ 1JZ = ^)3p, 



(3.3) 



where p = - [(ifa>) 4 g> (4) + (ifr>) 3 p t3) + (ifa)) 2 p (2) + {iw)p m + p (0) l- Expressions for coefficients p4,...,p t0) are 
presented in 13411 . 

Finally, time-dependent single-particle distribution function is obtained in the following form 



17™* =EC? op 



v T u 



H3^3t ^ fc=0 



fc„(fc) 



2 ' 4 

(ia») 5 + Z (iw) fc p fc 
fc=o 



Here C? are constant coefficients; t? are relaxation times represented by 

where q z are roots of the characteristics equation 

{q z f + p 4 tq z f + p 3 (q z f + p 2 {q z ) 2 + piiq z ) + p = 0. 



(3.4) 



(3.5) 



The dynamic dielectric susceptibility of a clamped crystal is defined as 



H m* „2 I(iw)V fc) „ 2 fl^[l(iw) fc p 
^ 3 KT)= lim 2 ^$£_ = ^ /J _« = i = L U^ 



fc„(fc) 



(iw) 5 + £ {ia>) k Pk 
k=o 



i=i 



The coefficients £3; are found from the following system of equations 



5 

! = 1 1 



13 i 



(3.6) 



f "11 


"12 


"13 


"14 


"15 


1 




( 131 ' 




' "1 I 


"21 


"22 


"23 


"24 


"25 






132 




"2 


"31 


"32 


"33 


"34 


"35 






133 




"3 


"41 


"42 


"43 


"44 


"45 






134 






< "51 


"52 


"53 


"54 


"55 


) 




I 135 j 




, "5 ; 



(3.7) 
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Here, the following notations are used 

"11 = T 2 T 3 T 4 T 5' "12 — T 1 T 3 T 4 T 5> "13 — T l T 2 T 4 T 5' w 14 = T 1 T 2 T 3 T 5' W 15 ~ T 1 T 2 T 3 T 4' 



"21 - T 2 T 3 T 4 +T 2 T 3 T 5 +T 2 T 4 T 5 +T 3 T 4 T 5' 
"23 = T 1 T 2 T 4 ~*~ T 1 T 2 T 5 T l T 4 T 5 T 2 T 4 T 5' 
"25 = T 1 T 2 T 3 T l T 2 T 4 T 1 T 3 T 4 T 2 T 3 T 4' 
"32 = tfjf + Tfrl + T^rf + t| + t| tf + t|t| ; 

H34 — T^T^ ^1^3 ^1^5 ^*2^3 ^2^5 ^*3^"5' 
n 41 = + T3 + T4 + Tg, n 42 = Tj + T3 + T4 + T5; 

«44 = Tj + + T3 + Tg, « 45 = Tj + T2 + T3 + T4; 



"22 = Tj T3 Tg + Tj T4T5 + Tj TjT| + Tf t|t|; 
"24 = Tf Tf T3 + Tf T§ X% + TjTjTf + Tf T3T5; 
"31 = TjTf + tf T4 + T3T4 + tf Tg + T3T5 + T^Tg J 
"33 = Tf T| + t\t\ + Tf Tf + x\ T* + T^T* + t£ T§ ; 

7^35 — "^2 1 ^3 ^ 1 ^"4 ^"2 ^3 ^2 ^4 ^3 ^"4 ' 

"43 =Tf + T2+T4 + T5, 

"51 = "52 = "53 = "54 = "55 = 1; 



"1 = £ £ n ^p (4) - "2 = f p n ^p (3) , " 3 = f /? n t*p 12) , 

i=l i=l i=l 

"4 = ft t?p (1) , " 5 = 4 T iP (0) - (3- 8 ) 

(=1 

The complex longitudinal dielectric permittivity of the deuteron subsystem of a mechanically 
clamped MD2XO4 crystal reads 

4(w, T) = l + 4^|' 3 Cw, T), e£(w, D = to^Cfti, T). 

A numerical analysis shows that the most important contribution to the dispersion of £33 (<y , T) is made by 
the first relaxational mode [X3(l) » %3(i)]. while the dispersion of the complex dielectric permittivity of a 
mechanically clamped crystal is close to the Debye one. If the piezoelectric coupling is omitted, £33 (w, T) 
transforms into the expression corresponding to Hill . 



4. Comparison of the numerical results with experimental data. 
Discussion 

Let us analyse the results of numerical calculations performed within the framework of the proposed 
model for longitudinal dynamic dielectric characteristics of the MCHi-^D^XC^ crystals and compare 
them with the corresponding experimental data. It should be noted that the theory developed in the 
previous sections, strictly speaking, is valid for the MD2XO4 type crystals only. The experimental data 
are available for the MCHi-jtD^XC^ crystals with different deuterations x (0 =S x *S 1). The experimen- 
tally established relaxational character of the dielectric dispersion of £33 (v, T) liiili^l in these crystals, 
as has been already mentioned, is associated with suppression of tunneling by short-range interactions. 
Therefore, we shall neglect the effects of proton tunneling in MCHi-^D^XC^. We shall assume that the 
proposed theory for these crystals is also valid if we use the averaged effective values of the model pa- 
rameters 

£(X) = £h(1 —X)+ £j)X, W{X) = Wh(1 ~ X) + WqX. 

In ji^], we calculated the static longitudinal, piezoelectric, elastic, and thermal characteristics of the 
MCHi-^D ^2X04 and explored their dependencies on the values of the model parameters. It was shown 
that at a proper choice of these values, a good quantitative agreement between the theoretical results and 
the corresponding experimental data was obtained. These sets of the model parameters are used herein 
in calculating the dynamic characteristics of MCHi-^D^XC^. 

The parameter a that sets the time scale of the dynamic processes in MtHi- ^D^XC^, is determined 
from the condition that theoretical results for frequency dependencies of £33 (v, T) at different tempera- 
tures agree with the experimental data. It is assumed that an weakly depends on temperature 

a= [p + r\at\)-w~ u , Ar=r-r c . 
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Table 1. The obtained optimum values of the model parameters for K(H 1 _ X D X )2P0 4 . 



X 


T c 
(K) 


To 
(K) 


£ 

fcB 

(K) 


w 

fcB 

(K) 


vaiO) 

fc'B 

(K) 


M3-,10" 18 
(esu-cm) 


M3 + ,10" 18 
(esu-cm) 


r° 

^33 


P- 
(s) 


R- 
(s/K) 


P + 
(s) 


R+ 
(s/K) 


0.00 


122.5 


122.5 


56.00 


422.0 


17.91 


1.46 


1.71 


0.73 


0.35 


0.0100 


0.43 


0.0160 


0.21 


146.0 


145.9 


63.78 


515.8 


23.18 


1.54 


1.79 


0.65 


0.85 


0.0095 


1.22 


0.0193 


0.29 


1SS.0 


154.8 


66.74 


551.5 


25.21 


1.57 


1.82 


0.62 


1.05 


0.0093 


1.51 


0.0217 


0.64 


191.0 


190.3 


79.71 


707.8 


32.34 


1.70 


1.96 


0.48 


1.76 


0.0385 


2.44 


0.0173 


0.79 


204.0 


203.1 


85.27 


774.8 


34.18 


1.76 


2.02 


0.42 


1.92 


0.0082 


2.65 


0.0151 


0.84 


208.0 


207.0 


87.12 


797.1 


34.63 


1.77 


2.03 


0.41 


2.02 


0.0081 


2.83 


0.0167 


0.91 


213.2 


212.2 


89.71 


828.4 


35.07 


1.80 


2.06 


0.38 


2.16 


0.0079 


2.88 


0.0130 


0.93 


215.0 


213.9 


90.45 


837.3 


35.36 


1.81 


2.07 


0.37 


2.20 


0.0079 


3.04 


0.0149 


0.99 


219.0 


217.9 


92.67 


864.1 


35.52 


1.83 


2.09 


0.35 


2.72 


0.0077 


4.21 


0.0189 


1.00 


220.1 


219.0 


93.05 


868.6 


35.76 


1.84 


2.10 


0.34 


2.84 


0.0077 


4.54 


0.0349 



X 


fcB 

(K) 


6 S 6 

fc B 

(K) 


t>n6 
fc B 

(K) 


OJ6 
fcB 

(K) 


4.10-10 
(dyn/cm 2 ) 


K 36 

(esu/cm 2 ) 


0.00 


-150.00 


82.00 


-500.00 


-400.00 


7.10 


1000.00 


0.64 


-142.73 


58.73 


-863.64 


-400.00 


6.59 


1727.27 


0.84 


-140.45 


51.45 


-977.27 


-400.00 


6.43 


1954.55 


0.93 


-139.43 


48.18 


-1028.41 


-400.00 


6.36 


2056.82 


1.00 


-138.64 


45.64 


-1068.18 


-400.00 


6.30 


2136.36 



Table 2. The obtained optimum values of the model parameters for RbH2P04 and KH2ASO4. 





T c 
(K) 


To 
(K) 


£ 

fc B 

(K) 


w 
fc B 
(K) 


fc'B 

(K) 


MS-JO" 18 
(esu-cm) 


M3 + ,10" 18 
(esu-cm) 


X33 


RbH 2 po 4 
KH 2 AS0 4 


147.6 
97.0 


147.6 
95.8 


60.00 
35.50 


440.0 
385.0 


29.13 
17.43 


1.50 
1.61 


2.00 
1.65 


0.40 
0.70 





Wa 
fc B 
(K) 


t>sB 
fcB 

(K) 


t>a8 
fcB 

(K) 


6'l6 
fcB 

(K) 


c 6 Vl0- 10 
(dyn/cm 2 ) 


K 36 

(esu/cm 2 ) 


RbH 2 P0 4 
KH 2 As0 4 


-130.00 
-170.00 


50.00 
130.00 


-500.00 
-500.00 


-300.00 
-500.00 


5.90 
7.50 


3000.00 
3000.00 





P- 


R- 


P + 


R + 




(s) 


(-) 


(s) 


(-) 


RM2PO4 


0.55 


0.0080 


0.93 


0.0140 


KH 2 As0 4 


0.47 


0.0160 


0.61 


0.0190 



The obtained optimum values of the model parameters are presented in table[T]for K(Hi_ x D x )2P0 4 
and in table[2]for RbH 2 P0 4 and KH 2 As0 4 . 

Note that fi3 + ,P+,R+ and /i3_,P_i?_ correspond to the paraelectric and ferroelectric phases, respec- 
tively. The temperature dependencies of the real and imaginary parts of the permittivity £33 (v, T) and 
£33 (v, T) at different frequencies for the KH2PO4, KD2PO4, RbH2P0 4 , and KH2ASO4 crystals are shown in 
figures|2}{5] Starting from a certain frequency v^, the low-frequency maximum in the temperature curve 
of £ 33 (v, T) is replaced with a sharp minimum at A T = K which widens and deepens with an increasing 
frequency that reaches £° 3 at v ~ 10 12 Hz. In KH 2 P0 4 v k = 33.2 GHz, in KD 2 P0 4 v k = 1.4 GHz, in RbH 2 P0 4 
v fc = 20.5 GHz, in KH 2 As0 4 v k = 20.8 GHz. The maximum of £ 33 (v, T) at A 7V, = \T n - T c \ decreases and 
smears out with an increasing frequency, whereas the magnitude of AT„ increases. With increasing fre- 
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-20 20 40 60 80 AT, K -20 20 40 60 80 AT, K 



Figure 2. The temperature dependence of e' 33 and e^, in KHyPC U at different frequencies v (GHz): 9.2 - 
1, (H; 33.2 - 2; 80 - 3; 154.2 - 4, □ 249 - 5, A (3J; 372 - 6, o\3^; 800 - 7. Symbols are experimental 
points; lines are the theoretical values. 




-20 20 40 60 80 AT, K -20 20 40 60 80 AT, K 



Figure 3. The temperature dependence of eL and e'^ in KD2PO4 at different frequencies v (GHz): 1.93 - 
1; 3.0 - 2, □ (jj; 10.0 -3,0 [37:]; 20.0 - 4, A [37|; 40.0 - 5; 80.0 - 6; 154.2 - 7. Symbols are experimental 
points; lines are the theoretical values. 
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-20 20 40 60 80 AT, K -20 20 40 60 80 AT, K 



Figure 4. The temperature dependence of eL and e'L in RbH2P0 4 at different frequencies v (GHz): 
0.25 - 1; 10.0 - 2; 27.0 - 3, V |i(i; 154.2 - 4, □ 250.2 - 5, A (ij; 372.0 - 6, o (ij; 700.0 - 7. Symbols are 
experimental points; lines are the theoretical values. 




Figure 5. The temperature dependence of e' 33 and in KH2ASO4 at different frequencies v (GHz): 9.2 - 
1, □ (H; 20.8 - 2; 80.0 - 3; 154.2 - 4, O («]]; 198.9 - 5, A 0; 250.2 - 6, > Q; 7, V gj; 372.0 - 8, < 0; 
700.0 - 9. Symbols are experimental points; lines are the theoretical values. 
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quency the magnitude of £33 ( v) decreases at all A T = \ T - T c \ . The maximal values of £ 33 (v) as well as the 
values of AT n are much larger in the paraelectric phase than in the ferroelectric phase. The dispersion 
of a real part of the permittivity £ 33 (v, T) in the ferrroelectric phase is observed in a narrow temperature 
range AT ~ 20 K, whereas in the paraelectric phase, AT is much larger, being of the order of 200 K. 

Let us note that taking into account the piezoelectric coupling, the calculated minimal values of £ 33 (v) 
at AT = at different frequencies are larger than those obtained within the model without the piezoelec- 
tric coupling. 

At a decreasing A T in the ferroelectric phase, the value of £ 3 ' 3 (v) increases, has a maximum at A T = 0, 
and decreases with an increasing AT in the paraelectric phase. At an increasing frequency, the maximal 
value of £33 (v) and the rate of its change with an increasing AT diminish. 

At v fc , the values of £ 33 (v, T) = £ 3 ' 3 (v, T) are 465 in KH 2 P0 4 , 520 in KD 2 P0 4 , 562 in RbH 2 P0 4 , and 330 
inKH 2 As0 4 . 

The proposed theory provides a good quantitative agreement with the experiment for KH 2 P0 4 (fig- 
ure[2) and a little worse agreement for the data of (37)] for KD 2 P0 4 (figure[3}, especially at AT < 20 K for 
£33 (v, T). However, it should, be noted that the values of £ 33 (v, T) obtained in (331 at frequencies above 
1 GHz have maxima at AT = K, rather than minima. 

The temperature dependence of £33 (v, T) in RbH 2 P0 4 measured in |41] is appropriately and well de- 
scribed by the present theory, except for the values of £ 33 (v, T) at v = 154.2 GHz and A T < 20 K (figure[4}. 
The theory and experimental data of [40] for £33 (v, T) at v = 27 GHz are also in a good agreement. The 
obtained theoretical results for £33 (v, T) at 198 and 366 are only in qualitative agreement with the data 
of |{426l , which, in their turn, are in disagreement with the results of other measurements of j4lh . 

The calculated temperature dependencies of £33 (v, T) and £33 (v, T) accord well with the ones mea- 
sured in (4lll for KH 2 As0 4 at different frequencies starting from the submillimeter range (figure[5}. The 
data for £33 (v, T) obtained in (36fl at v = 9.2 GHz are in a somewhat worse agreement with the theory, 
especially at AT < 20 K. 

Figures [6lfl3l contain the calculated temperature dependencies of the real and imaginary parts of 
longitudinal dynamic dielectric permittivity £ 33 (v, T) and £ 3 ' 3 (v, T) of clamped K(Hi_ x D x ) 2 P0 4 crystals at 
different deuterations x and frequencies along with the corresponding experimental data. 




60 
50 
40 
30 
20 
10 



E 33 




V 




□\ 















20 40 60 80 100120 AT, K 





-40-20 
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Figure 6. The temperature dependence of e 33 and £33 in K(Ho.7gDo. 2 i)2P0 4 at different frequencies 
v (GHz) 154.2 - 1, □; 249.0 - 2, A; 372.0 - 3, o. Symbols are experimental points; lines are the 
theoretical values. 

With an increasing deuteration x in K(Hi_j;D x ) 2 P0 4 , the magnitude of £ 33 (v, T) decreases, whereas 
AT n increases. 

At an isomorphic replacement K — ► Rb, P— ► As, the maximal values of £ 33 (v, T) remain almost un- 
changed, whereas AT n slightly increase. 

It should be noted that in the MH 2 X0 4 , the experimental data of (H, HH HH correspond to the 
region of dielectric permittivity dispersion. At the same time, for KD 2 P0 4 in the measurements of (iill . 
the submillimeter frequencies correspond to a high-frequency "tail" of the dispersion, whereas in the 
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Figure 7. The temperature dependence of e' 33 and in K(Ho.36Do.64)2 po 4 at different frequencies 
v (GHz) I35J: 154.2 - 1, □; 249.0 - 2, A; 372.0 - 3, o. Symbols are experimental points; lines are the 
theoretical values. 




-40 -20 20 40 60 80 100 120 AT, K -40 -20 20 40 60 80 100 120 AT, K 



Figure 8. The temperature dependence of e' 33 and in K(Ho.i6Do.84)2P04 at different frequencies 
v (GHz) 3.">|: 154.2 - 1, □; 249.0 - 2, A; 372.0 - 3, o. Symbols are experimental points; lines are the 
theoretical values. 




-40 -20 20 40 60 80 100 120 AT, K -40 -20 20 40 60 80 100 120 AT, K 



Figure 9. The temperature dependence of e' 33 and in K(Hrj.o7Do.93)2P04 at different frequencies 
v (GHz) (Hi: 154.2 - 1, □; 249.0 - 2, A; 372.0 - 3, o. Symbols are experimental points; lines are the 
theoretical values. 



33705-13 



R.R. Levitskii, I.R. Zachek, A.S. Vdovych 




Figure 10. The temperature dependence of £ 33 and e 3 ' 3 in K(Hi_ x D x )2P04 at v=9.2 GHz and for 
different x lial : 0.0 - 1, o; 0.29 - 2, •; 0.99 - 3, •. Symbols are experimental points; lines are the theo- 
retical values. 




-40 -20 20 40 60 80 100 120 AT, K -40 -20 20 40 60 80 100 120 AT, K 



Figure 11. The temperature dependence of e 33 and e 33 in KOHi-j^D^PC^ at v=138.6 GHz and for differ- 
ent x |43j]: 0.63 - 1, O; 0.91 - 2, □. Symbols are experimental points; lines are the theoretical values. 




Figure 12. The temperature dependence of e 33 and e 33 in K(Ho.22Do.78)2 po 4 at different frequencies 
v (GHz) |38]: 8.6 - 1, A; 9.7 - 2, □; 26.5 - 3, o. Symbols are experimental points; lines are the theoretical 
values. 



data of I37I1 . this is the low-frequency tail. Further experimental measurements of £33 (v, T) at v > 10 are 
required to evaluate the validity of the calculated [v, T) . 

The most graphic illustration of the dispersion of the real and imaginary parts of the dielectric permit- 
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tivity £33(11*, T) in MCHi-^D^XC^ would be their frequency-temperature plots drawn in wide frequency 
and temperature ranges. Such plots for theoretical dependencies along with the experimental points are 
presented in figures Il4l Il5l for K(Ho.rj7D .93) 2 P0 4 , in figures H6lH7l for RbH 2 P0 4> and in figures IT81H91 for 
KH 2 As0 4 . 




are experimental points; lines are the theoretical values. 



Let us analyse the changes in the real and imaginary parts of £33 (<y) in the MCHi-^D^XC^ crystals 
at replacing H -* D, K -* Rb, and P —> As. At AT - +0 K, the dispersion frequency [i.e., the frequency 
of the maximum of e^M] is 33.2 in KH 2 P0 4 , 1.93 in KD 2 P0 4 , 20.5 in RbH 2 P0 4 , and 20.8 in KH 2 As0 4 . 
The linewidth [i.e., the difference between frequencies of the maximum and half maximum of £ 3 ' 3 (<y)] is 
12.0 GHz in KH 2 P0 4 , 4.8 in KD 2 P0 4 , 7.8 in RbH 2 P0 4 , and 7.2 in KH 2 As0 4 . 

The temperature dependencies of the inverse relaxation time (t^) -1 in K(Hi_ x D A -) 2 P0 4 along with the 
values estimated from different experimental measurements are presented in figure [20] The calculated 
values of the relaxation times t^ 3 4 , in contrast to Tj, are unlikely to depend on temperature and are 
much smaller than the value of Tj. The theory provides a satisfactory agreement with the experiment for 
temperature curves of the relaxation time. A certain difference between the relaxation times estimated 
from the dielectric permittivity and ultrasound measurements is due to the contributions into attenuation 
from the mechanisms irrelevant for the permittivity (e.g., scattering by admixtures). 
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Figure 15. The frequency-temperature dependence of £ 3 ' 3 in K(Ho.o7Do.93)2P°4- ■> A > • - I3al . Symbols 
are experimental points; lines are the theoretical values. 




Figure 16. The frequency-temperature dependence of £33 in RbH2P04. - |4(jl; A, o - (4l||. Symbols are 
experimental points; lines are the theoretical values. 




Figure 17. The frequency-temperature dependence of £33 in RbH2P04. 0- (43]; □, a, o- (ij. Symbols are 
experimental points; lines are the theoretical values. 



33705-16 



Longitudinal relaxation of mechanically clamped crystals 




Figure 19. The frequency-temperature dependence of £33 in KH2ASO4. □ - (i^]; o, a, >, v, < - lilll . 
Symbols are experimental points; lines are the theoretical values. 




120 140 160 180 200 220 240 260 T, K 



Figure 20. The temperature dependence of the inverse polarization relaxation time at different x: 0.0 - 
1, > jH, <J3, O 0.07 - 2, □ del; 0.21 - 3, A \s&; 0.43 - 4, V \s&; 0.72 - 5, ■ (1(1; 0.805 - 6, ▲ (H; 
0.84 - 7, O l52ll : 0.93 - 8, ▼ (i^l; 1.0 - 9, • H3|. Symbols are experimental points; lines are the theoretical 
values. 
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5. Conclusions 

Within the framework of the modified proton ordering model, taking into account a linear over the 
strain £q contribution into the energy of the proton subsystem, and using the four-particle cluster approx- 
imation, we calculate the longitudinal dynamic characteristics of mechanically clamped crystals of the 
KH2PO4 family. For the partially deuterated crystals M(Hi_ x D x )2X04, these characteristics are obtained 
within the mean crystal approximation. The data for £33 (v, T) and c 3 ' 3 (v, T) presented by different groups 
of experimentalists are analyzed and systematized. At the proper choice of the theory parameters for the 
MCHi-^D.^XC^ crystals, we obtain a good quantitative description of the available experimental data 
for £33 (v, T) and c 3 ' 3 (y, T) . For the first time, the dispersion of the longitudinal dynamic dielectric permit- 
tivity of clamped crystals of the KH2PO4 family is explored in wide temperature and frequency ranges. It 
should be noted that the effect of piezoelectric coupling on the dielectric characteristics of these crystals 
is essential. In the present paper, the observed temperature behavior of £ 33 (v, T) in the phase transition 
region at different frequencies has been appropriately described for the first time. 
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l"lo3flOB>KHfl pe/iaKcai^in MexamMHO 3aTi/icHyTi/ix Kpi/iCTa/iiB Ti/my 
KH 2 P0 4 

P.P. /leBi/iLjbKi/iiP, LP. 3aHeP, A.C. BflOBi/i ^ 

iHCTHTyT (f>i3HKW KOHfleHCOBamix cucTeM HAH YKpaiHn syn. CB£HL4ii4bKoro, 1, JlbBiB, 79011, YKpaiHa, 
Hai^ioHa^bHww yHiBepcuTeT "ZlbBiBCbKa no/iiTexHiKa" syn. C. BaHflepw 12, 79013, ZlbBiB, YKpaiHa 

y paMKax MOflwcfiiKOBaHoT MOfle/ii npoTOHHOro BnopHflKyBaHHH cer'HeToaKTWBHWx KpncTa/iiB ciM'l KH2PO4 3 Bpa- 
xyBaHH^M ^iHitiHoro 3a fle^opMai^ieio BHecKy £6 b eHepriio npoTOHHOi' cucTeMw b Ha6/in>KeHHi 40Tnpn4acTHH- 
koboto K/iacrepa b Mexax flWHaMHHOi MOfle/ii nnay6epa OTpnMaHO Bnpa3 fl/ia no3flOBXHboT flWHaMHHOi flie/ie- 
KTpn4H0T npoHUKHOCTi MexaHiHHO 3a™cHyToro Kpucra/iy. npw Ha^e>KHOMy Bn6opi napaMeTpiB Teopi'i OTpwMa- 
ho flo6pnw Ki/ibKicHufi onnc HaflBHnx eKcnepwMeHTa/ibHwx flaHnx fl/ifl u,vtx Kpuda/iiB. 

K/iHJMOBi c/ioBa: cemeToeneKTpi/iKvi, K/iacrepHe Ha6fli/i)KeHHn, flie/ieKTpn^Ha npoHUKHicTb, naciA penaKcau,i'i 
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